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Abstract 

We give an explicit algorithm and source code for combining alpha streams 
via bounded regression. In practical applications typically there is insufficient 
history to compute a sample covariance matrix (SCM) for a large number of 
alphas. To compute alpha allocation weights, one then resorts to (weighted) 
regression over SCM principal components. Regression often produces alpha 
weights with insufficient diversification and/or skewed distribution against, 
e.g., turnover. This can be rectified by imposing bounds on alpha weights 
within the regression procedure. Bounded regression can also be applied to 
stock and other asset portfolio construction. We discuss illustrative examples. 
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1 Introduction 


With technological advances there is an ever increasing number of alpha streams]^ 
Many of these alphas are ephemeral, with relatively short lifespans. As a result, 
in practical applications typically there is insufficient history to compute a sample 
covariance matrix (SCM) for a large number of alpha streams - SCM is singular. 
Therefore, directly using SCM in, say, alpha portfolio optimization is not a optional 

One approach to circumvent this difficulty is to build a factor model for alpha 
streams [56]. Because of the utmost secrecy in the alpha business, such factor models 
must be build in-house - there are no commercial providers of “standardized” factor 
models for alpha streams. As with factor models for equities, such model building 
for alphas requires certain nontrivial expertise and time expenditure. 

Therefore, in practice, one often takes a simpler path. As was discussed in 
more detail in [56] . one can deform SCM such that it is nonsingular, and then use 
the so-deformed SCM in, say, Sharpe ratio optimization for a portfolio of alphas. 
For small deformations this then reduces to a cross-sectional weighted regression 
of the alpha stream expected returns [55] . The regression weights are the inverse 
sample variances of the alphas. The columns of the loadings matrix, over which the 
expected returns are regressed, are nothing but the first K principal components of 
SCM corresponding to its positive (be., non-vanishing) eigenvalues [ 56] . 

Regression often produces alpha weights with insufficient diversification and/or 
skewed distribution against, e.g., turnover. Thus, if some expected returns are 
skewed, then, despite nonunit regression weights (which suppress more volatile al¬ 
phas), the corresponding alpha weights can be larger than desired by diversification 
considerations. Also, the principal components know nothing about quantities such 
as turnover I/] A simple way of obtaining a more “well-rounded” portfolio composi¬ 
tion is to set bounds on alpha weights. This is the approach we discuss here. 

When individual alpha streams are traded on separate execution platforms, the 
alpha weights are non-negative. By combining and trading multiple alpha streams 
on the same execution platform - the framework we adapt here - one saves on 
transaction costs by internally crossing trades between different alpha streams (as 
opposed to going to the market) 0 Then the alpha weights can be negative. 

When alpha weights can take both positive and negative values, the bounded 
regression problem simplifies. It boils down to an iterative algorithm we discuss in 
Section [2] This algorithm can actually be derived from an optimization algorithm 


3 Here “alpha” - following the common trader lingo - generally means any reasonable “expected 
return” that one may wish to trade on and is not necessarily the same as the “academic” alpha. 
In practice, often the detailed information about how alphas are constructed may not be available, 
e.g., the only data available could be the position data, so “alpha” then is a set of instructions to 
achieve certain stock holdings by some times t\, t 2 > ■ • • 

4 For a partial list of hedge fund literature, see, e.g., ra-Eni and references therein. For a partial 
list of portfolio optimization and related literature, see, e.g., EU-Ea and references therein. 

5 One approach to rectify this is to add a turnover-based factor to the loadings matrix [SB] . 

6 For a recent discussion, see [57] . 
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with bounds (in a factor model context) discussed in [58] by taking the regression 
limit of optimization. We also give R source code for the bounded regression algo¬ 
rithm in Appendix [A] Appendix [B] contains some lcgalese. We conclude in Section 
[31 where we also discuss bounded regression with transaction costs following [59] . 


2 Bounded Regression 

2.1 Notations 

We have N alphas a*, i = Each alpha is actually a time series a?j(t s ), 

s — 0,1,..., M, where t 0 is the most recent time. Below a* refers to ai(t 0 ). 

Let Cij be the sample covariance matrix (SCM) of the N time series ai(t s ). If 
M < N, then only M eigenvalues of C\j are non-zero, while the remainder have 
“small” values, which are zeros distorted by computational rounding^ 

Alphas ctj are combined with weights Wi. Any leverage is included in the defini¬ 
tion of he., if a given alpha labeled by j G {1,..., N} before leverage is a' (this 
is a raw, unlevered alpha) and the corresponding leverage is Lj : 1, then we define 
aj = Lj a'j. With this definition, the weights satisfy the condition 

N 

XlH = 1 (!) 

i =1 

Here we allow the weights to be negative as we are interested in the case where the 
alphas are traded on the same execution platform and trades between alphas are 
crossed, so one is actually trading the combined alpha a = YliLi a* wy. 

2.2 Weighted Regression 

When SCM Cij is singular and no other matrix ( e.g ., a factor model) to replace it 
is available, one can deform SCM such that it is nonsingular, and then use the so- 
deformed SCM in, say, Sharpe ratio optimization for a portfolio of alphas [55]. For 
small deformations this reduces to a cross-sectional weighted regression of the alpha 
stream expected returns [56] • The regression weights z t (not to be confused with 
the alpha weights W\) are the inverse sample variances of the alphas: Zi = 1/Cu. 
The columns of the loadings matrix Aj^, A = 1,..., K, over which the expected 
returns are regressed, are nothing but the first K principal components of SCM 
corresponding to its positive (he., non-vanishing) eigenvalues. However, for now we 
will keep A % a general (e.g., one may wish to include other risk factors in A^ [56]). 

' Actually, this assumes that there are no N/As in any of the alpha time series. If some or all 
alpha time series contain N/As in non-uniform manner and the correlation matrix is computed by 
omitting such pair-wise N/As, then the resulting correlation matrix may have negative eigenvalues 
that are not zeros distorted by computational rounding. 
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The weights uy are given by: 


Wi = l Zi £i ( 2 ) 

where £ t are the residuals of the cross-sectional regression of ay over (without 
the intercept, unless the intercept is subsumed in A^, that is - see below) with the 
regression weights zf 


N K 

£% Qy ^ ^ Zj Oij ^ ) A-iA A. jb Q AB (3) 

j =1 a,b =1 


where Qjf B is the inverse of 


N 

QaB = Zi A iA hiB (4) 

i =1 

and the overall factor 7 in (J2]) is fixed via (JT|). Note that we have 

N 

VA G {1, ..., K} : ^2w i A iA = 0 (5) 

i=l 

So, the weights Wi are neutral w.r.t. the risk factors defined by the columns of the 
loadings matrix A.^. 


2.3 Bounds 


Since the weights wy can have either sign, we will assume that the lower and upper 
bounds on the weights 


satisfy the conditions 


< Wi < wf 

(6) 

wi < 0 

(T) 

wf > 0 

(8) 

wi < wf 

(9) 


The last condition is not restrictive: if for some alpha labeled by i we have w ~ = 
w 4 + , then we can simply set wy = w~ and altogether exclude this alpha from the 
bounded regression procedure below. Also, if, for whatever reason, we wish to have 
no upper/lower bound for a given uy, we can simply set wf = ± 1 . 

The bounds can be imposed for diversification purposes: e.g., one may wish to 
require that no alpha has a weight greater than some fixed (small) percentile £, i.e., 
\wi\ < £, so wf = ±£. One may also wish to suppress the contributions of high 
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turnover alphas, e.g., by requiring that |wj| < £ if t* > r*, where ly is the turnover J^] 
r* is some cut-off turnover, and £ is some (small) percentile. Bounds can also be 
used to limit the weights of low capacity^ alphas. Etcfi^ 

2.4 Running a Bounded Regression 

So, how do we impose the bounds in the context of a regression? There are two 
subtleties here. First, we wish to preserve the factor neutrality property (]5]h which 
is invariant under the simultaneous rescalings w, t —> (wi (where ( is a constant). If 
we simply set some uy to their upper or lower bounds, this generally will ruin the 
rescaling invariance, so the property ([5]) will be lost. Second, we must preserve the 
normalization condition ([1]). In fact, it is precisely this normalization condition that 
allows to meaningfully set the bounds w±, as the regression itself does not fix the 
overall normalization coefficient 7 in ([ 2 ]) , owing to the rescaling invariance wy —?• (yuy. 

Here we discuss the bounded regression algorithm. To save space, we skip the 
detailed derivation as it follows straightforwardly by taking the regression limit 
of optimization with bounds in the context of a factor model, both of which are 
discussed in detail in [58] 1^1 

Let us define the following subsets of the index i e J = {1,..., N}: 


Wi = wf, i E J + (10) 

Wi = w~, i £ J~ (11) 

7 = J + U J- (12) 

I=J\J (13) 

Further, let 

on = 7 a i (14) 

y A = ^Zi ai A iA + ^ w t w i A '^ ( 15 ) 

ieJ ieJ+ i&J- 


8 Here the turnover (over a given period, e.g., daily turnover) is defined as the ratio t 7 ; = Di/Ii 
of total dollars Di (long plus short) traded by the alpha labeled by i over the corresponding total 
dollar holdings L (long plus short). 

9 By capacity I* for a given alpha we mean the value of the investment level L for which the 
P&L Pi(Ii) is maximized (considering nonlinear effects of impact). 

10 Since the regression we consider here is weighted with the regression weights Zi = l/Cu, this 
already controls exposure to alpha volatility, so imposing bounds based on volatility would make 
a difference only if one wishes to further suppress volatile alphas. 

11 The regression limit of optimization essentially amounts to the limit £ z 2 = y £ z 2 , 77 —» 0, 
£ 2 = fixed, where is the specific (idiosyncratic) risk in the factor model with the factor loadings 
matrix identified with the regression loadings matrix (and the K x K factor covariance matrix 
becomes immaterial in the regression limit) - see [58| for details. 
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where 7 is to be determined (see below). Then we have 


K 


Wi = Zi ai 


E A*a Qab Vb ) , i E J 


a,b =1 


K 


Vi E J + : Zi 


OLi 


E Qab ) — 


M); 


A,B=1 

K 


Vi E J : Zj | a* — E Aja Qab 2/b ) < 


- 


A,B=1 


( 16 ) 

(17) 

(18) 


where 0 1 is the inverse of the K x K matrix 0: 

Qab = ^ A*a Ajb (19) 

ieJ 

Here the loadings matrix A*a must be such that Q is invertiblef^l Also, note that 
Wi, i E J given by (fT6jl together with Wi = wf, i E J + and w t = w~, i E J~ satisfy 
dHJ), as they should. 

Note that, for a given value of 7 , (fl5|) solves for i/a given J + and J~. On the 
other hand, (1171) and (fT 8 l) determine J + and J~ in terms of i/a- The entire system is 
then solved iteratively, where at the initial iteration one takes J ® = J, so that J +< -°) 
and are empty. However, we still need to fix 7 . This is done via a separate 

iterative procedure, which we describe below. 

Because we have two iterations, to guarantee (rapid) convergence, the J ± it¬ 
eration (that is, for a given value of 7 ) can be done as follows. Let wf s) be such 
that 


Vi E J : w t < w \ s ■* < wf 


N 

VAe{1,...,K}: E A*a = 0 

i= 1 


( 20 ) 

( 21 ) 


At the (s + l)-th iteration, let wf +l> be given by (TT 6 |) for i E j( s \ with wf +1) = wf 
for i E . This solution satisfies (J5D, but may not satisfy the bounds. Let 


- „>+!) 


'(s) 


hi(t) = wf'* +t q i: t 6 [ 0 , 1 ] 


( 22 ) 

(23) 


Then 

wf +1) = hi(U) = w\ s) + U qi (24) 

12 This is the case if the columns of A*a are comprised of the first I\ principal components of 
SCM Cij corresponding to its positive eigenvalues. However, as mentioned above, here we keep 
the loadings matrix general. 
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where i* is the maximal value of t such that hi(t) satishes the bounds. We have: 


Qi > 0 : pi = min 
Qi < 0 : pi = max 


(«+i) 

5 

(»+i) 


£* = min 


Pi ~ w. 
qi 


(«) 


7^ 0) f G t/ 


Now, at each step, instead of (H7li and (fT 8 jh we can define J ± h+ 1 ) v j a 

Vi 6 J+<‘ +1 > : = w+ 

Vi 6 J- (,+11 : fii, ( ” +1) = w~ 


(25) 

(26) 

(27) 


(28) 

(29) 


where y) ' +l ' is computed iteratively as above and we can take wP = 0 at the initial 
iteration. Unlike era and (EEBD , (|28l) and ()29l) add new elements to J ± one (or a 
few) element (s) at each iteration. 

The convergence criteria are given by 


j+(s+i) _ j+{s ) 

j~( s + 1) — j-(s) 


(30) 

(31) 


These criteria are based on discrete quantities and are unaffected by computational 
(machine) precision effects. However, in practice the equalities in (1251) and (129|) 
are understood within some tolerance (or machine precision) - see the R code in 
Appendix [A] We will denote the value of u)j s+1 ' 1 at the final iteration (for a given 
value of 7 , that is) via w t . 

Finally, 7 is determined via another iterative procedure as follows (we use su¬ 
perscript a for the 7 iterations to distinguish it from the superscript s for the ,J ± 
iterations): 


^( a+1 ) 


ry(. a ) 


z^i =1 



(32) 


where w^ is computed as above for 7 = 7 ^. To achieve rapid convergence, the 
initial value 7 *-°) can be set as follows: 


7 


(o) = 



(33) 


where e % are the residuals of the weighted regression (without bounds) given by (ED • 
The convergence criterion for the 7 iteration is given by 


(o+l) _ (a) 


( 34 ) 


understood within some preset computational tolerance (or machine precision). 
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The R code for the above algorithm with some additional explanatory documen¬ 
tation is given in Appendix [Aj Note that this code is not written to be “fancy” or 
optimized for speed or in any other way. Instead, its sole purpose is to illustrate 
the bounded regression algorithm as it is described above in a simple-to-understand 
fashion. Some legalese relating to this code is given in Appendix |B] 

2.5 Application to Stock Portfolios 

Above we discussed the bounded regression algorithm in the context of computing 
weights for portfolios of alpha streams. However, the algorithm is quite general and 
- with appropriate notational identifications - can be applied to portfolios of stocks 
or other suitable instruments. In fact, it can also be applied outside of finance. 
Here, for the sake of definiteness, we will focus on stock portfolios, in fact, we will 
assume that they are dollar neutral, so both long and short positions are a 1 lowed [^1 

2.5.1 Establishing Trades 

Let us first discuss establishing trades, i.e., we start from nil positions and establish 
a portfolio of N stocks. Instead of alpha streams, our index i e {1,..., N} = J now 
labels the stocks. We will denote the desired dollar (not share) holdings via Hi, and 
the total dollar investment (long plus short) via I: 

N 

/ = E m ( 35 ) 

i =1 

Let Wi = Hi/I. These are now our stock weights (analogous to the alpha weights). 
Then we have the familiar normalization condition 

N 

y \w/\ = 1 (36) 

1=1 

However, normally, one imposes bounds on Hi, not on Wi. For example, in the case 
of establishing trades one may wish to cap the positions such that: i) not more 
than a small percentile £ of the total dollar investment / is allocated to any given 
stock - this is a diversification constraint; and ii) only a small percentile £ of ADDV 
(average daily dollar volume) V t is traded - this is a liquidity constraint (see below). 
In this case we have the following bounds on the dollar holdings H,;. 

H~ <H t <Hf (37) 

Hf = ±min ({ /, f V() (38) 

In this case the upper and lower bounds are symmetrical. In some cases, such as 
for hard-to-borrow-stocks, we may have some H~ = 0. In other cases one may not 

13 Various generalizations are possible, some more straightforward than others. 
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wish to have a long position in some stocks. Etc. We will only assume that H~ < 0 
and > 0, in line with our discussion above for the bounds on the weights, which 
are then given by 

wf = Ht/1 ( 39 ) 

The final touch then is that instead of cq one uses some expected returns Ei in the 
case of stocks. The rest goes through exactly as above for a suitably chosen A^. 

2.5.2 Rebalancing Trades 

With rebalancing trades, we have the current dollar holdings H* and the desired 
dollar holdings Hi. In this case, one may wish to cap the positions such that: i) not 
more than a small percentile £ of the total dollar investment / is allocated to any 
given stock - this is the same diversification constraint as above; ii) only a small 
percentile £ of ADDV Vi is traded - this the same liquidity constraint as above; and 
iii) not more than a small percentile £' of ADDV V is allocated to any given stock 
- this is another liquidity constraint stemming from the consideration that, if the 
portfolio must be liquidated swiftly (e.g., due to an unforeseen event), to mitigate 
liquidation costs, the positions are capped based on liquidity. Here £' typically can 
be several times larger than £ - the portfolio can be built up in stages as long at 
each stage the bounds are satisfied. The bounds on Hi now read: 

\H t \ < mm (£/, i'Vi) (40) 

\H t - H*\ < £ V; (41) 

It is more convenient to rewrite these bounds in terms of the traded dollar amounts 
Di = Hi - H*-. 


D~ < Di < Dt (42) 

Dt = min (min (£ /, £' Vi) - H*, £ V-) > 0 (43) 

D~ = max (-min (£ /, £' ly) - //*. - £ V:) < 0 (44) 

and we are assuming that \H*\ < min (£ /, £' V). Furthermore, we will assume that 
H* itself satisfies (0: 

N 

VAg{1,...,A'}: H i = 0 (45) 

i =1 


Then the bounded regression algorithm can be straightforwardly applied to the 
weights Wi and aq defined as follows: 


Wi = Hi/1 
Xi = Di/1 


(46) 

(47) 


In the iteration we now use Xi instead of Wi, while in the 7 iteration we still use w t . 
Then the rest of the algorithm goes through unchanged. Let us note, however, that 
the source code given in Appendix [X] is written with alpha weights in mind, so while 
it can be adapted to the case of stock portfolios in the case of establishing trades, 
straightforward modifications are required to accommodate rebalancing trades. 

2.5.3 Examples: Intraday Mean-Reversion Alphas 

To illustrate the use of the algorithm, we have employed it to construct portfolios 
for intraday mean-reversion alphas with the loadings matrix in the following 5 
incarnations: i) intercept only (so K = 1); ii) BICS (Bloomberg Industry Classifi¬ 
cation System) sectors; iii) BICS industries; iv) BICS sub-industries; and v) the 4 
style factors prc, mom, hlv and vol of [60j plus BICS sub-industries. The regres¬ 
sion weights are the inverse sample variances: z* = 1 / C tl (see below). I 11 the cases 
ii)-v) above the intercept is subsumed in the loadings matrix Aj^. Indeed, we have 
AjA = 1, where G is the set of columns of A t A corresponding to sectors in 
the case ii), industries in the case iii), and sub-industries in the cases iv) and v). 
Consequently, the resultant portfolios are automatically dollar neutral. 

The portfolio construction and backtesting are identical to those in |61j . where 
more detailed discussion can be found, so to save space, here we will only give a brief 
summary. The portfolios are assumed to be established at the open and liquidated 
at the close on the same day, so they are purely intraday and the algorithm of Section 
12.5.11 for the establishing trades applies. The expected returns A) for each date are 
taken to be E l = —A,, where Ri = In (p° pen j and for each date p° pen i s 
today’s open, while Pf ose is yesterday’s close adjusted for splits and dividends if the 
ex-date is today. So, these are intraday mean-reversion alphas. 

The universe is top 2000 by ADDV V), where ADDV is computed based on 21- 
trading-day rolling periods. However, the universe is not rebalanced daily, but also 
every 21 trading days (see [61] for details). The sample variances C %% are computed 
based on the same 21 -trading-day rolling periods, and are not applied daily, but 
also every 21 trading days, same as the universe rebalancing (see [ 6 TJ for details). 
We run our simulations over a period of 5 years (more precisely, 252 x 5 trading 
days going back from 9/5/2014, inclusive). The annualized return-on-capital (ROC) 
is computed as average daily P&L divided by the total (long plus short) intraday 
investment level / (with no leverage) and multiplied by 252. The annualized Sharpe 
Ratio (SR) is computed as the daily Sharpe ratio multiplied by \/252. Cents-per- 
share (CPS) is computed as the total P&L divided by the total shares traded. On 
each day the total (establishing plus liquidating) shares traded for each stock are 
given by Qi = 2\H i \/P° pen (see [61] for details). 

For comparison purposes, the results for regressions without bounds are given in 
Table [H The results for the bounded regressions, with the bounds on the desired 
holdings set as 

\H t \ < 0.01 K (48) 
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so not more than 1% of each stock’s ADDV is bought or sold, are given in Table [2] 
and the corresponding P&Ls are plotted in Figure 1. Thus, as expected, adding the 
liquidity bounds has the diversification effect on the portfolios, so the Sharpe ratios 
are substantially improved - as usual, at the expense of (slightly) lowering paper 
ROC and CPS. Note that, even with tight liquidity bounds, the 4 style factors prc, 
mom, hlv and vol of [60J add value, further validating the 4-factor model of |60j . 

3 Concluding Remarks 

One - but not the only - way to think about bounded regression is as an alternative 
to optimization with bounds when the latter is not attainable. In fact, as mentioned 
above, bounded regression is a zero specific risk limit of optimization with bounds in 
the context of a factor model. So, when a factor model is not available, e.g., in the 
context of alpha streams, bounded regression can be used in lieu of optimization. 

In this regard, one can further augment the bounded regression algorithm we 
discussed above by including linear transaction costs, as in [59], A systematic ap¬ 
proach is to start with optimization with bounds and linear transaction costs in the 
context of a factor model as in [5E] and take a zero specific risk limit. Non-linear 
transaction costs (impact) in the context of alpha weights can be treated using the 
approximation discussed in [59] using the spectral model of turnover reduction [El]. 


A The R Code 

Below we give R (R Package for Statistical Computing, http://www.r-project.org) 
source code for the bounded regression algorithm we discuss in the main text. 
The entry function is calc.bounded.lm(), which runs the 7 iteration loop and 
calls the function bounded.lm() , which runs the J ± iteration loop. The args() 
of calc.bounded.lm() are: ret, which is the iV-vector of alphas cp (or, more gener¬ 
ally, some other returns); load, which is the N x K loadings matrix A iAl weights, 
which is the iV-vector of the regression weights zf, upper, which is the iV-vector of 
the upper bounds w +: lower, which is the N -vector of the lower bounds w”; and 
prec, which is the desired precision with which the output weights the iV-vector 
of which calc.bounded.lm() returns, must satisfy the normalization condition (P). 
Internally, bounded.lm() calls the function calc.bounds (), which computes wj s+1 
in (|24|) at each iteration. The code is straightforwardly self-explanatory. Jp, Jm 
in bounded.lm() correspond to J ± . One subtlety is that, when restricting Aj^ to 
J C J, in the case of binary industry classification (e.g., when A*^ corresponds 
to BIOS sub-industries, which can be small), the so-restricted A iA may have null 
columns, which must be omitted and the code below does just that. For uou-biuary 
cases, one may wish to augment the code to ensure that the matrix Q <- t (load [Jt, 
take]) w.loadfJt, take] is nonsingular (and if it is, then remove the culprit 
columns in A iA or otherwise modify the latter); however, for non-binary A iA and 
generic regression weights this should not occur except for special, non-generic cases. 
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calc.bounded.lm <- function(ret, load, weights, upper, lower, prec = le-5) 

{ 

reg <- lm(ret ~ -1 + load, weights = weights) 
x <- weights * residuals(reg) 
ret <- ret / sum(abs(x)) 

repeatj 

x <- bounded.lm(ret, load, weights, upper, lower) 
if(abs(sum(abs(x)) - 1) < prec) 
break 

ret <- ret / sum(abs(x)) 

} 


return(x) 

} 


bounded.lm <- function(ret, load, 

{ 

calc.bounds <- function(z, x) 

{ 


} 


q <- x - z 

p <- rep(NA, length(x)) 
pp <- pmin(x, upper) 
pm <- pmax(x, lower) 
p[q > 0 ] <- pp[q > 0 ] 
p [q < 0] <- pm[q < 0 ] 
t <- (p - z)/q 
t <- min(t, na.rm = T) 
z <- z + t * q 
return(z) 


weights, 


upper, lower, tol = le-6) 


if(!is.matrix(load)) 

load <- matrix(load, length(load), 1) 

n <- nrow(load) 
k <- ncol(load) 

ret <- matrix(ret, n, 1) 
upper <- matrix(upper, n, 1) 
lower <- matrix(lower, n, 1) 
z <- diag(weights) 


11 



w.load <- z %*% load 
w.ret <- z %*% ret 

J <- rep(T, n) 

Jp <- rep(F, n) 

Jm <- rep(F, n) 
z <- rep(0, n) 

repeatj 

Jt <- J & !Jp & !Jm 
y <- t(w.load[Jt, ]) 7o*7» ret[Jt, ] 
if(sum(Jp) > 1) 

y <- y + t (load [Jp, ]) 7»*7» upper [Jp, ] 
else if(sum(Jp) == 1) 

y <- y + upper[Jp, ] * matrix(load[Jp, ], k, 1) 
if(sum(Jm) > 1) 

y <- y + t (load [Jin, ]) 7o*7o lower [Jm, ] 
else if(sum(Jm) == 1) 

y <- y + lower[Jm, ] * matrix(load[Jm, ], k, 1) 
if(k > 1) 

take <- colSums(abs(load[Jt, ])) >0 
else 

take <- T 

Q <- t(load[Jt, take]) 7»*7o w.load[Jt, take] 

Q <- solve(Q) 
v <- Q 7»*7o y[take] 

xJp <- Jp 
xJm <- Jm 

x <- w.ret - w.load[, take] 7 o*7o v 
x[Jp, ] <- upper [Jp, ] 
x[Jm, ] <- lower [Jm, ] 

z <- calc.bounds(z, x) 

Jp <- abs(z - upper) < tol 
Jm <- abs(z - lower) < tol 

if(all(Jp == xJp) & all(Jm == xJm)) 
break 

} 


return(z) 


12 



} 


B DISCLAIMERS 

Wherever the context so requires, the masculine gender includes the feminine and/or 
neuter, and the singular form includes the plural and vice versa. The author of this 
paper (“Author”) and his affiliates including without limitation Quantigic® Solu¬ 
tions LLC (“Author’s Affiliates” or “his Affiliates”) make no implied or express 
warranties or any other representations whatsoever, including without limitation 
implied warranties of merchantability and fitness for a particular purpose, in con¬ 
nection with or with regard to the content of this paper including without limitation 
any code or algorithms contained herein (“Content”). 

The reader may use the Content solely at his/her/its own risk and the reader 
shall have no claims whatsoever against the Author or his Affiliates and the Author 
and his Affiliates shall have no liability whatsoever to the reader or any third party 
whatsoever for any loss, expense, opportunity cost, damages or any other adverse 
effects whatsoever relating to or arising from the use of the Content by the reader 
including without any limitation whatsoever: any direct, indirect, incidental, spe¬ 
cial, consequential or any other damages incurred by the reader, however caused 
and under any theory of liability; any loss of profit (whether incurred directly or 
indirectly), any loss of goodwill or reputation, any loss of data suffered, cost of pro¬ 
curement of substitute goods or services, or any other tangible or intangible loss; 
any reliance placed by the reader on the completeness, accuracy or existence of the 
Content or any other effect of using the Content; and any and all other adversities 
or negative effects the reader might encounter in using the Content irrespective of 
whether the Author or his Affiliates is or are or should have been aware of such 
adversities or negative effects. 

The R code included in Appendix [A] hereof is part of the copyrighted R code 
of Quantigic® Solutions LLC and is provided herein with the express permission of 
Quantigic® Solutions LLC. The copyright owner retains all rights, title and interest 
in and to its copyrighted source code included in Appendix [A] hereof and any and 
all copyrights therefor. 
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Table 1: Simulation results for the 5 alphas via regression without bounds discussed 
in Section 12.5.31 


Alpha 

ROC 

SR 

CPS 

Regression: Intercept only 

33.59% 

5.59 

1.38 

Regression: BICS Sectors 

39.28% 

7.05 

1.61 

Regression: BICS Industries 

42.66% 

8.19 

1.75 

Regression: BICS Sub-industries 

45.25% 

9.22 

1.84 

Regression: 4 Style Factors plus BICS Sub-industries 

46.60% 

9.85 

1.90 


Table 2: Simulation results for the 5 alphas via bounded regression discussed in 
Section 12.5.31 


Alpha 

ROC 

SR 

CPS 

Regression: Intercept only 

29.66% 

7.36 

1.25 

Regression: BICS Sectors 

35.32% 

9.89 

1.48 

Regression: BICS Industries 

39.25% 

12.00 

1.65 

Regression: BICS Sub-industries 

42.23% 

14.13 

1.75 

Regression: 4 Style Factors plus BICS Sub-industries 

43.70% 

15.54 

1.82 
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Figure 1. P&L graphs for the intraday alphas discussed in Section T2.5.3L with a summary 
in Table [2j Bottom-to-top-performing: i) Regression over intercept only, ii) regression over 
BIOS sectors, iii) regression over BICS industries, iv) regression over BICS sub-industries, 
and v) regression over 4 style factors prc, mom, hlv and vol of |60| plus BICS sub-industries. 
The investment level is $10M long plus $10M short. 
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